Inhomogeneous phase formation on the border of itinerant ferromagnetism 
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A variety of analytical techniques suggest that quantum fluctuations lead to a fundamental insta- 
bility of the Fermi liquid that drives ferromagnetic transitions flrst order at low temperatures. We 
present both analytical and numerical evidence that, driven by the same quantum fluctuations, this 
first order transition is pre-empted by the formation of an inhomogeneous magnetic phase. This 
occurs in a manner that is closely analogous to the formation of the inhomogeneous superconduct- 
ing Fulde-Ferrel-Larkin-Ovchinnikov state. We derive these results from a field theoretical approach 
supplemented with numerical Quantum Monte Carlo simulations. Our analytical approach repre- 
sents a considerable simplification over diagrammatic methods, makes contact with older analyses 
of the unitarity limit, and enables a simple physical picture to emerge. 
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Many magnetic materials display second order ferro- 
magnetic phase transitions. The transition temperature 

[ can be tuned using external parameters such as doping 
and pressure. Hertz realized that tuning such a transi- 
tion to zero temperature could give rise to a new type of 
critical phenomena, for which he coined the term quan- 
tum criticality This led to a tremendous experimen- 

] tal and theoretical effort that has had some notable suc- 
cesses with the universal scaling predicted for the quan- 

\ tum critical regime being seen in a handful of materi- 
als Q. However, in the majority of systems currently in- 
vestigated, new behavior intervenes before the transition 
temperature can be tuned to zero. In many cases, such 

• as ZrZna 3, UGc2 [3, MnSi [|, and C0S2 the second 
order transition becomes first order before the quantum 
critical point is reached. Moreover, recent experimental 

\ evidence points to phenomena that go beyond the first 
order transition; materials such as ZrZn2 UGe2 01 j 

; Ca3Ru207 H, NbFe2 H, and Sr3Ru207 aU display 
unusual behavior in the vicinity of the putative quantum 
critical point. 

This failure to find a naked quantum critical point has 
lead to speculation that it represents a fundamental prin- 
ciple [Tlj . Diagrammatic calculations that extend beyond 
the standard Moriya-Hertz-Millis [l| theory of itinerant 
quantum criticality suggest a breakdown of the Landau 
expansion around the quantum critical point (T^ . [isl . [l^ . 
This raises the question of how to connect the diagram- 
matic calculations to the well-established second order 
perturbation approach [l^. This approach accounts for 
all orders of vacuum scattering amplitude, and predicts 
that the itinerant ferromagnetic transition is first order 
at low temperature. 

We show that, when a linearization of the electron dis- 
persion about the Fermi surface is permissible, the first 
order magnetic transition is pre-empted by the forma- 
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tion of an inhomogeneous magnetic phase in a manner 
closely analogous to the inhomogeneous superconducting 
Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state [l^- Our 
results are consistent with the effects of non-analyticities 
that appear in extensions to the Hertz-Millis theory [I^, 
[T3| . Here, we adopt an alternative field theoretical ap- 
proach [l7| to provide analytical evidence for our pic- 
ture. As well as resolving the connection between the 
second order perturbation theory approach and the di- 
agrammatic analysis, it provides insight into how quan- 
tum fluctuations drive the reconstruction of the phase 
diagram. Our analytical considerations are supported by 
Quantum Monte Carlo (QMC) simulations. 

This letter is divided into two main parts: We first 
develop our analytical results and then discuss the QMC 
simulations that support them. Since one of the main 
advantages of our approach is the simplification that it 
affords, we outline the main steps of our analysis in their 
entirety. We follow this with a discussion of the effect 
of quantum fluctuations upon the phase diagram of the 
homogeneous fcrromagnet - showing both how older re- 
sults may be recovered from our analysis and the non- 
analyticities revealed by diagrammatics. Next, we show 
how the same fluctuations can drive a spatial modula- 
tion. After presenting our QMC results, we conclude 
with a discussion of possible extensions to the work. 

In metals, the long-range component of the Coulomb 
interaction is screened. We therefore take as our start- 
ing point a free electron system interacting through 
a contact (Hubbard-like) repulsive interaction, g [11]. 
The corresponding partition function may be expressed 
as a fcrmionic coherent state path integral Z = 
rpj. g-/3(H-MW) ^ J x)^,e-s with the action 

IY.^'- (^r + C) j 5^+V^-^-V^+ • (1) 

Here / = /g^ dr / Sx, Cfc = ^ where Cj^ denotes the 
dispersion and /i represents the chemical potential. To 
develop an effective Landau theory of the magnetic tran- 
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sition, Hertz introduced a scalar Hubbard-Stratonovich 
decoupling of the interaction in the spin channel 
However, this form of decoupling neglects the potential 
impact of soft transverse field fluctuations. It is the 
latter that are responsible for driving the second order 
transition flrst order and in turn promote the instabil- 
ity towards inhomogcneous phase formation. Therefore, 
we will introduce a general Hubbard-Stratonovich decou- 
pling which incorporates fluctuations in all of the spin $ 
and charge p sectors. 

Defining # = m + 0, where m denotes the puta- 
tive saddle-point value of the magnetization field and 
4> the fluctuations, one obtains the partition function 
Z = J V\l)'DpV4>c~^ where the action now takes the form 

S ^ S g[m^+4>''-p')+i^[dr+C+gp-g(T-{rn+(i>)]i: 0. 
We wish to study the potential for spatially modulated 
order where the magnetization forms a conical spiral, ro- 
tating about some axis (set along z) with a pitch vector 
q. To simplify the analysis, it is helpful to transform to a 



rotating basis in which the magnetization becomes spa- 
tially uniform, m = (to±, 0, Setting ip i— > e^'^'^'^^/'^ip^ 

\ / cos(q • r) - sin(q • r) \ ( <l^^\ 
(j>y \ ^ \ sin(q-r) cos(q • r) 0y , 

0z / V 1 y V 0, y 

and integrating over the fcrmion degrees of freedom yields 
the action 

S = j g{TT? + 02 _ p2) 

- Tr \vL\dr + Ck+a,q/2 - ffCT • m +gp - gCT • 0] . (2) 

V ' 

G-1 

Then, integrating over p and an expansion of the ac- 
tion to second order in g leads to 



Z, = cxp 



grn^ + TrlnG"^ -Tr.gn+" - -Tr52(n+-n-+ - n++fl"-) 



where H''* = G^G" denotes the Lindhard function and = [dr -I- — u) ^ with 

k.q ' 



ek-hq/2 + ek-q/2 _^ V(ek-q/2 - ek-q/2 + '^9'^^)'^ + (25?71||)2 



(3) 



r 



representing the energy of the electrons in plane-wave 
states with momentum k and spin-up or down relative 
to the mean-field spiral. To remove the unphysical ul- 
traviolet divergence due to the contact nature of the 
interaction potential and arising from the term in the 
action second order in g, we must affect the standard 
regularization of the linear term TrgiH^ setting g 

2/cFa/^z.-2(2/cFa/7r^.)Vl/EL3,4««,q)'^(^k4,q)«,q + 



e^^ q ^ ej|^^ q — e^^ q)^^ UM^ where the prime indicates 
that the summation is subject to the momentum conser- 
vation ki -|- k2 = -|- k4. This regularization allows us 
to characterize the strength of the interaction through 
the dimensionless parameter kpa, where kp denotes the 
Fermi wave vector, a is the s-wave scattering length, and 
1/ is the density of states at the Fermi surface. 

Finally, carrying out the Matsubara summations, one 
obtains the following expression for the free energy: 



k,s 



k,q" l^k q; 



^N+N- 




f de+de 



- £+)/?„ (-k,e ) 



"(4i,q)"(ek2,q) 



ki,q ' '^k2,q k3,q "^k4,q 



— , (4) 



where n(e) = 1/(1-1- e'^'^'^ ''^) is the Fermi distribution, 
^q = Ek'^(4,q)' and p±(k,e) = Ek' "(e^'+k/2.q)[l - 

"(«^'-k/2,q)]^(« - ^k'+k/2,q + 4-k/2,q) 



1/(1 + e'3(^ 

q = Ek'^(4,q)' and 

^ ' ' the spm- 

up(down) particle-hole density of states. 

Spatially Uniform Case: To orient our discussion, we 
first consider the implications of the fiuctuation correc- 
tions on the magnetic phase diagram for a free particle 
dispersion ek = /2m* without accounting for spatial 
modulation. In this case Eq. ([4]) reduces to that ob- 
tained in Ref. [lB| from second order perturbation the- 
ory. From this result, one finds that fiuctuations drive the 
second order ferromagnetic transition flrst order at tem- 
peratures below that of the tricritical point. Tt ~ 0.2Tf, 
where Tp is the Fermi temperature (see Fig. [T}. Sub- 
stituting the low energy form of the particle-hole den- 
sity of states, Pq^Q(k, e) = ed{kkp — k'^/2 — e)/27rfc, 

where fcp = •y/2m*(// ± 2kpam/TTi/), into the fiuctuation 
correction to the free energy, and expanding in powers 
of magnetization, one recovers a singular term of order 
In arising from particle-hole excitations with mo- 
mentum near to 2fcF, the same non-analyticity as was 
found diagrammatically in Refs. [l^, The forma- 

tion of a finite magnetization increases the phase-space 
available for the formation of virtual intermediate pairs 
of particle- hole pairs, and this phase space enhancement 
ultimately drives the ferromagnetic transition flrst order. 
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Spatially Modulated Case: Quantum fluctuations also 
lead to spatial modulation of the magnetic order and 
further reconstruction of the magnetic phase diagram. 
The results of this analysis are shown in Fig. [1] which 
is obtained from Eq. ([4]) through a Landau expansion 
in ni [20|. At low temperatures and kpa = 0.84, a sec- 
ond order transition into an inhomogencous spin phase 
with pitch O.lfcp pre-empts the first order phase tran- 
sition at fcpa = 1.055. The expansion cannot however 
describe how far the textured phase penetrates into the 
uniform ferromagnetic phase. For a quadratic disper- 
sion, the magnetization and spiral wave-vector enter the 
single electron energy in the combination (k • q/m*)^ + 
{4:kFamj_/Tn')'^ , where we have restricted to a planar spi- 
ral since it always has lower energy than a conical spi- 
ral at zero magnetic field. Upon linearizing the electron 
energy at the Fermi surface, and to^ enter the free 
energy in the same way (up to the angular factors that 
accompany q and which are integrated over). Therefore, 
the spatial modulation enters as if it where a direction 
dependent magnetization. As a consequence, the coeffi- 
cient of in our Taylor expansion is proportional to 
the coefficient of q^m^. When the m^-term becomes 
negative - and a first order transition becomes favorable 
- the q^TO^-term also becomes negative favoring a spa- 
tial modulation that emerges from the tricritical point. 
This isprecisely the same situation as seen in the FFLO 
state 

To assess the validity of the perturbative scheme, we 
turn now to the numerical Quantum Monte Carlo anal- 
ysis of the Stoner Hamiltonian ([1]) making use of the 
CASINO program [23|- These methods are based upon 
optimizing a trial wave function and are restricted to zero 
temperature. Our approach mirrors that used in previous 
studies of itinerant ferromagnetism [23,[H,[23|. The vari- 
ational wave function used in our simulation, ip = Dc""^ , 
is a product of a Slater determinant, D, that takes ac- 
count of the Fermion statistics and occupation of single 
particle orbitals, and a Jastrow factor, J, that accounts 
for electron correlations. 

The Slater determinant consists of plane-wave spinor 
orbitals containing both spin-up and spin-down electrons, 
D — det({V'kGfcT , V'kefe; })■ Although not an eigcnstatc of 
total spin, the constituent spin states of least weight pro- 
vide the dominant contribution to the variational state 
energy. In the case of uniform magnetization, for com- 
putational efficiency, we factorize the Slater determinant 
into an up and a down-spin determinant [2l|. The spin 
textured phase is described by non-coUincar spins, which 
have only recently been studied within the Variational 
Monte Carlo (VMC) These studies lead us to de- 

scribe a planar spin spiral with a trial wave function 
that contains the spinors ipk = e"i''"/^(e''' '', e"* "") and 
-(/ik = e~"' '"/^(— e''' '', e^'*' '") which explicitly fix the spin 
spiral orientation. For q = this would recover the 
factorized form for the Slater determinant employed in 
the uniform case. Our simulations are carried out in a 
unit cell with periodic boundary conditions commensu- 
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FIG. 1: The lower phase diagram shows the first order (dashed 
line) and second order (upper solid line) ferromagnetic order- 
ing if the system were restricted to be homogeneous. The 
first order transition is pre-empted by a second order phase 
transition into a modulated ferromagnetic state bounded by 
the solid line. The dotted line shows where the uniform ferro- 
magnetic phase would form, if the transition were restricted 
to be second order; this putative phase boundary not only 
comes after the first order transition and joins smoothly to 
the second order phase boundary at temperatures above the 
tricritical point, but also has negative slope at T = 0, consis- 
tent with Refs. The upper graph shows the wave vector 
q of the modulated phase boundary. 



rate with the pitch of the spiral. 

The Jastrow factor, J, accounts for electron-electron 
correlations. It consists of polynomial and plane-wave 
expansions in electron-electron separation, and a poly- 
nomial backflow function. In the spiral case, the Jastrow 
factor is restricted to be spin independent to maintain the 
spin spiral orientation and the wave function antisymme- 
try. The trial wave functions were optimized using QMC 
methods. In the uniform case, the optimization was per- 
formed in two steps using VMC and Diffusion Monte 
Carlo, whereas only VMC calculations were performed 
for the textured state. To model the repulsive contact 
potential between the electrons we employ the modified 
Poschl- Teller interaction [2^ which has smooth edges so 
that the QMC configurations can sample it faithfully [13 . 

Firstly, constraining the magnetization to be spatially 
uniform, an estimate of the ground state magnetization 
for different interaction strengths recovers the expected 
phase diagram (Fig. [2]), revealing a first order phase tran- 
sition into the itinerant magnetic phase. This provides a 
platform upon which to construct the full textured phase 
diagram. Allowing for spatial modulation of the magne- 
tization, we find that an inhomogencous magnetic phase 
pre-empts the transition into the uniform phase. The re- 
sulting textured phase has similar extent and wave vector 
to the analytical prediction lending support to the con- 
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FIG. 2; The lower panel shows the variation of the ground 
state magnetization M with interaction strength fcpa at zero 
temperature. The dashed line corresponds to the uniform 
phase, and the solid line is the textured phase. The upper 
panel shows the wave vector of the inhomogeneous magnetic 
phase, and the discrete values of q sampled in the investigation 
are highlighted by the horizontal dotted lines. 



elusions of the perturbative field theoretic analysis. 

In conclusion, quantum fluctuations are known to drive 
the itinerant ferromagnetic transition first order. We 
have shown that the same mechanisms lead to the devel- 



opment of inhomogeneous magnetic order in the vicinity 
of the tricritical point, resulting in a phase diagram which 
mirrors that of the superconducting FFLO phase. Our 
results are consistent with recent diagrammatic analyses 
and reveal the connection between these works ITsI 
and older second order perturbation calculations [l5|. 
Moreover, our approach represents a considerable analyt- 
ical simplification and more clearly reveals the underlying 
physical processes. 

There are several directions in which this analysis 
might be extended. An electronic band dispersion can 
drive a similar reconstruction of the electronic phase di- 
agram [2^ . It would be informative to investigate the 
interplay between such effects and quantum fluctuations. 
As well as the possibility of spatially modulated mag- 
netism, recent works have suggested that a band dis- 
persion might lead to d-wave distortion of the Fermi 
surface (29j . Such electron nematics may be viewed as 
melted versions of the spatially modulated magnetism 
that we consider. Indeed, quantum fluctuations might 
drive the formation of an electron nematic phase in itin- 
erant ferromagnets. 
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